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The two-photon exchange contribution to the single spin asymmetries with the spin orientation 
normal to the reaction plane is discussed for elastic electron-proton scattering in Regge regime. 
For this, the equivalent photon approximation is adopted. In this case, hadronic part of the two- 
photon exchange amplitude describes real Compton scattering (RCS). The imaginary part of helicity 
conserving RCS amplitudes are related to total photoabsorption cross section. The contribution of 
the photon helicity nipping amplitudes is estimated by the two pion exchange in the i-channel. We 
observe the double logarithmic enhancement in the electron mass but find it's contribution to the 
asymmetry negligibly small in the forward kinematics. These results are in strong disagreement 
with the existing calculation. 
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I. INTRODUCTION 

Recently, the new polarization transfer data for the 
electromagnetic form factors ratio Ge/Gm Q raised a 
lot of interest for the two photon exchange (TPE) physics 
in elastic electron proton scattering. These new data ap- 
peared to be incompatible with the Rosenbluth data 
A possible way to reconcile the two data sets was pro- 
posed |3j, which consists in a more precise account on 
the TPE amplitude, the real part of which enters the ra- 
diative corrections to the cross section (Rosenbluth) and 
the polarization cross section ratio in a different manner. 
At present, only the IR divergent part of the two pho- 
ton exchange contribution, corresponding to one of the 
exchanged photon beeing soft, is taken into the experi- 
mental analysis £| . Two model calculations exist for the 
real part of the TPE amplitude pj, 0| , and they quali- 
tatively confirm that the exchange of two hard photons 
may be responsible for this discrepancy. In order to ex- 
tract the electric form factor in the model independent 
way, one has thus to study the general case of Compton 
scattering with two spacelike photons. These two photon 
contributions are important for the electroweak sector, as 
well. 

In view of this interest, parity-conserving single spin 
asymmetries in elastic ep-scattering with the spin orien- 
tation normal to the reaction plane regain an attention. 
These observables are directly related to the imaginary 
part of the TPE amplitude and have been an object of 
theoretical studies in the 1960's and 70's 0. By ana- 
lyticity, the real part of the TPE amplitude is given by 
a dispersion integral over its imaginary part. Therefore, 
a good understanding of this class of observables is ab- 
solutely necessary. Recently, first measurements of the 
beam normal spin asymmetry Bn 1 have been performed 
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1 In the literature, also the A n notation for beam normal spin 
asymmetry or vector analyzing power was adopted. 



in different kinematics 

Though small (several to tens ppm), this asymmetry 
can be measured with a precision of fractions of ppm. 
Before implementing different models for the real part, 
where an additional uncertainty comes from the disper- 
sion integral over the imaginary part, one should check 
the level of understanding of the imaginary part of TPE. 
These checks have been done for the existing data. Inclu- 
sion of the elastic (nucleon) intermediate state only j^j led 
to negative asymmetry of several ppm in the kinematics 
of SAMPLE experiment but was not enough to describe 
the data. The description of the beam normal spin asym- 
metry within a phenomcnological model which uses the 
full set of the single pion electroproduction ^1 did not 
give satisfactory description at any of the available kine- 
matics. Especially intriguing appears the situation with 
the SAMPLE data with electron lab energy Ei ab = 200 
MeV, which is just above the pion production threshold 
where the theoretical input is well understood. On the 
other hand, an EFT calculation without dynamical pions 
was somewhat surprisingly very successful in describ- 
ing this kinematics for B n . This success suggests that to 
the given order of chiral perturbation theory, the role of 
the dynamical pions for this observable might be not too 
large. 

Though even at low energies the situation with the 
imaginary part of TPE amplitude is by far not clear, 
an attention has to be paid to high energies, as well, 
since the dispersion integral which would give us its real 
part, should be performed over the full energy range. 
Due to relative ease of measuring B n within the frame- 
work of parity violating electron scattering, new data 
from running and upcoming experiments [12| will stim- 
ulate further theoretical investigations of this new ob- 
servable. A calculation of B n in hard kinematics regime 
at high energy and momentum transfer was performed 
recently in the framework of generalized parton distribu- 
tions (GPD's) and resulted in asymmetries of ~ 1.5 ppm. 

Since a ppm effect measurement at high momentum 
transfers is an extremely complicated task, we concen- 
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trate in this work on the forward kinematics. For this 
kinematics, a calculation exists 0], where an observa- 
tion is made that the contribution of the situation where 
the exchanged photons are nearly real and overtake the 
external electron kinematics is enhanced as ln 2 (— t/m 2 ), 
with m the electron mass and t < the clastic momen- 
tum transfer. The authors of take the forward VCS 
amplitude and DIS structure functions as an input, and 
provide an estimation of B n in this kinematics as large as 
25-35 ppm. We demonstrate that such a treatment is not 
adequate and leads to an overestimation of the effect by 
an order of magnitude for the JLab energy Ei a b — 5.77 
GeV. In order to provide an estimate of the contribu- 
tion of this kinematics to B n , we use the most general 
real Compton scattering amplitude. We observe that the 
double logarithmic enhancement is not a consequence of 
taking the forward limit but a pure kinematical effect 
within the loop integral. We demonstrate that this treat- 
ment of the problem leads to much smaller values for this 
asymmetry than those quoted in [T^. We provide fur- 
thermore a clear explanation of this discrepancy. In fact, 
it is a consequence of the well known fact that the forward 
doubly virtual Compton scattering amplitude has a non- 
analytic behaviour in the real photon point Q 2 — 0, with 
Q 2 the virtuality of the incoming and outgoing photons. 
By now, it has been noticed for the low energy limit of 
the real and (doubly) virtual Compton scattering, where 
even the leading terms in the low energy expansion are 
sensitive to the order of taking either the low energy limit 
v — > or Q 2 — > first [ljj. Our observation generalizes 
this non-analytical behaviour from the kinematical point 
v = 0, t = to the whole forward axis t = 0. 

The article is organized as follows: in Sectional we de- 
fine the kinematics, general ep-scattering amplitude and 
the observables of interest; in Section ITTT1 the two pho- 
ton exchange mechanism and the photons kinematics is 
studied; the equivalent photons or quasi real Compton 
scattering approximation and its implementation for the 
case of B n is given in Section IIVI we present our results 
in Section which are followed up by a discussion and a 
summary. 



II. ELASTIC ep-SCATTERING AMPLITUDE 



In this work, we consider elastic electron-proton scat- 
tering process e(k) + p(p) — > e(fc') + p(p') for which we 
define: 



and choose the invariants t — q 2 < 2 and v = (P ■ K) /M 
as the independent variables. M denotes the nucleon 
mass. They are related to the Mandelstam variables s = 
(p + k) 2 and u = (p — k') 2 through s — u = Mvlv and 
s + u + t = 2M 2 . For convenience, we also introduce 
the usual polarization parameter e of the virtual photon, 
which can be related to the invariants v and t (beglecting 
the electron mass m): 



v 2 - M 2 t{1 + r) 



(2) 



v 2 + M 2 t(1 +t)' 

with r = — t/(4M 2 ). Elastic scattering of two spin 
1/2 particles is described by six independent amplitudes. 
Three of them do not flip the electron helicity [3j, 

e 2 _ 

T no fu P = —u(k )^fj,u(k) (3) 



u(p') ( Gm7 



F 2 — 

M 



3 M 2 



while the other three are electron helicity flipping and 
thus have in general the order of the electron mass m 
M: 



Tfup - — — 



m e 
M-t 



u{k')u{k) ■ u{p') [Fi 
{- F 6 u(k')-f 5 u(k) ■ u(p'h 5 u(p) 



f4)u( P ) 



(4) 



In the one-photon exchange (Born) approximation, two 
of the six amplitudes match with the electromagnetic 
form factors, 

cfr>,i) = c M {t), 

Ff°™(M) = F 2 (t), 

FiZ, n Mt) = (5) 

where Gm (t) and F 2 (t) are the magnetic and Pauli form 
factors, respectively. For further convenience we define 
also Ge = Gm — (1 + i~)F 2 and Fx = Gm ~ F 2 which 
in the Born approximation reduce to electric form factor 
Ge and Dirac form factor F\, respectively. For a beam 
polarized normal to the scattering plane, one can define 
a single spin asymmetry, 



B,. 



<7f + <T L 



(6) 



where cr-f (crjj denotes the cross sesction for an unpolar- 
ized target and for an electron beam spin parallel (an- 
tiparallel) to the normal polarization vector defined as 



sz = o, 



[k x k'} 
\kxk'\ 



(7) 



P = 
K = 



P + P 
2 

k + k' 



k — k! = p' — p, 



(1) 



In elastic ep-scattering, the usual notation for the momentum 



transfer is Q 



-q 2 but we prefer the more general notation t 



to avoid confusion with the incoming and outgoing photon vir- 
tualities in forward doubly virtual Compton scattering we will 
be concerned in the following. 



3 



normalized to (S ■ S) — — 1. Similarly, one defines the 
target normal spin asymmetry A n . It has been shown 
m the early 1970's that such asymmetries are directly 
related to the imaginary part of the T-matrix. Since 
the electromagnetic form factors and the one-photon ex- 
change amplitude are purely real, B n obtains its finite 
contribution to leading order in the electromagnetic con- 
stant a em from an interference between the Born ampli- 
tude and the imaginary part of the two-photon exchange 
amplitude. In terms of the amplitudes of Eas. (|3l4|) . the 
beam normal spin asymmetry is given by: 



B n = ~ — V^O^Tk/IT^ (tG 



eGi 



tG m ImF 3 



G E InLF 4 + F^lmF^Y (8) 



and using the cm. frame of the electron and proton, one 
has: 



with \k\ = 



Ql,2 

s-M 2 
2V^ 



= 2|fc||fci|(l -COS 01,2), (11) 

k the three momentum of the incom- 



-w 2 -\-rn 2 \ 2 



) 2 — m 2 



ing (and outgoing) eletron, \ki\ — y y — 

that of the intermediate electron, and cos 82 = 
cos 6 cos 0i + sin sin 0i cos 0. The kinematically al- 
lowed values of the virtualities of the exchanged photons 
(the restriction is due to the fact that the intermediate 
electron is on-shell) are represented by the internal area 
of the ellypses shown in Fig. [21 The ellypses are drawn 

Accessible values of Q, 2 , Q, 2 



For completeness, we also give here the expression of 
target normal spin asymmetry T„ 3 in terms of invariant 
amplitudes: 



T n = y/2e{l + e)^ (tG 2 m + eG\) 1 



(9) 



+ 



F1ITD.F2 — F2Im.F1 
v 



^-Ge-G m \-I^F, 



III. TWO PHOTON EXCHANGE 




FIG. 1: Two-photon exchange diagram. 

The imaginary part of the two-photon exchange (TPE) 
graph in Figgis given by 



ImAl27 = 



^d|fcikjflfa 
2£i(27r) 3 



QlQ\ 



:W^(w 2 ,Qi,> 



11), 



(10) 



where W fl ' / (w 2 , Q 2 , Q 2 ,) is the imaginary part of doubly 
virtual Compton scattering tensor. Q 2 and Q 2 denote 
the virtualities of the exchanged photons in the TPE di- 
agram, and w is the invariant mass of the intermediate 
hadronic system. We next study the kinematics of the 
exchanged photons. Neglecting the small electron mass 




FIG. 2: Kinematically allowed values of the photon vitualities 



inside a square whose side is defined through the external 
kinematics (k) and the invariant mass of the intermedi- 
ate hadronic state (w 2 or ki), while the form solely by 
the scattering angle. If choosing higher values of the 
mass of the hadronic system w 2 < s, it leads to scal- 
ing the size of the ellypse by a factor of s s _^ 2 . In the 
limit w 2 — (y/s — m e ) 2 , the ellypses shrink to a point 
at the origin and both photons are nearly real. This is 
not a collinear singularity, however, since the real pho- 
tons' energy remains large enough in order to provide the 
transition from nucleon with the mass M to the interme- 
diate state X with the mass w. Instead, the intermediate 
electron is soft, fcf ~ (m e , 0), therefore this kind of kine- 
matics does not lead to an IR divergency which can only 
occur if the intermediate hadronic state is the nucleon 
itself. In the following we are going to study this kine- 
matical situation in more detail. 



IV. QUASI-RCS APPROXIMATION 



3 Also A n notation for target normal spin asymmetry exists in the 
literature. 



We consider the kinematical factors under the integral 
over the electron phase space of Ea. (|lU|) at the upper 
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limit of the integration over the invariant mass of the 
intermediate hadronic state, w — > \fs — m e : 



1 



k\ 1 

Ex Q\Q\ ~ 4k 2 E, 



1 1 

m 4fc 2 



(12) 



In this range of the hadronic mass w, the exchanged pho- 
tons are real, and the contribution of real Compton scat- 
tering (RCS) to the 27-exchange graph is enhanced by a 
factor of 1/m (it is not a singularity since B n has a factor 
of m in front). This enhancement, however, only appears 
for the beam asymmetry B n since the target asymmetry 
does not involve the electron spin flip and the electron 
mass can be neglected in the electron propagator fa + to 
in the numerator of Ea. (|10|l . The remaining fa cancels 
this 1/m behaviour. 

We next rewrite the hadronic tensor in Eq. I|10[l identi- 
cally as 

W^{W 2 ,Ql,Q 2 2 ) = W^{ Sl 0,0) (13) 

+ (W^(W 2 ,QlQ 2 2 )-W^(s,0,0)), 

so that only the first term is now enhanced in this limit, 
while the second term vanishes at the RCS point by 
construction. Equivalent photon or quasi-real Compton 
scattering (QRCS) approximation consists in assuming 
the first term to be dominant due to the kinematical 
enhancement under the integral and in neglecting the 
second one. The question of the validity of such an ap- 
proximation should be raised for the beam normal spin 
asymmetry. We will come back to this point in Section 
IVl However, in general, this kind of contributions com- 
ing from QRCS kinematics will allways be present in the 
full calculation, since the second term in Ea. (|14|) is con- 
structed in such a way that the resulting integral is reg- 
ular at the QRCS point. In the following, the QRCS 
approximation will be used. Hence, the hadronic ten- 
sor can be taken out of the integration over the electron 
phase space. The remaining integrals are 



d 3 fa 



1 



2£i(27r) 3 Q\Ql 
d 3 fa k{ 



We next list the result of the integration in the limit 
of high energies (i.e., neglecting terms ~ m n /E): 
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hp 

IlK 
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s- M 2 (4E 2 
16tt 2 M 4 - su 11 \~t 



(15) 



1 



-4tt 2 < 



In 



IE 



AMv 
+ -—hp 



T° 



4Mv 



t 



tip- 



For the details, we address the reader to the Appendix. 
The three integrals Iq, I® k , and hp obviously classify as 



log 2 to, log 1 to, and log to, respectively. In Ref. |14| only 
the first integral was calculated. 

The RCS tensor may be taken for instance in the basis 
of Prange 0] or, equivalently of Berg and Lindner [l7j . 



RCS 



( p'M p'v 



if J "n v 



-(S3 + fBi 



P a n 2 



P' 2 n 2 



N, 



(16) 



with the vectors defined as P' = P — ^j^-K, = 
E^ vot ^P v K a q & such that P' ■ K = P' ■ n = n ■ K = 0. 
The amplitudes Bi are functions of v and t. This form of 
Compton tensor is convenient due to the simple form 
of the tensors appearing in Eq. ((TE|) . However, the 
amplitudes Bi possess kinematical singularities and con- 
straints, therefore in the following we will also use an- 
other set of invariant amplitudes introduced as combi- 
nations of BiS in This latter set of amplitudes is 
widely used in the dispersion analysis of Compton exper- 
iments, high energy contributions to them are well stud- 
ied, and we will take an advantage of this knowledge. 

Before contracting the leptonic and hadronic tensors, 
we notice that the amplitude B-j can only contribute to 
the invariant amplitude Fq, since it contains the N'j^N 
structure. F$ does not contribute at leading order in 
to to neither obserevable of interest, therefore B7 will 
be neglected in the following. The remaining tensors in 
Eq. (|16|) are symmetric in indices [iv. 



lmM^ RCS 



= e 2 u'^ v (^i + TOJoh^Wncs 
= u^-fa + mlotyW^sgpu 
+ 2W R ^ cs hf 1 u-^ v u. 



(17) 



Finally, we identify different terms in Eq. 1)17(1 with the 
structures, together with amplitudes Fi v ..6 parametriz- 
ing elastic ep-scattering amplitude in Eqs. (USUI) and after 
some algebra we find for the invariant amplitudes for the 
elastic electron-proton scattering in the QRCS approxi- 
mation: 



Iihf: 



QRCS 



-2th P ImB 6 
-MtI 1P 



Im 



Bi — B 3 



2Mt 



M 4 



ImK 



QRCS _ 



-MHhplm 



Im 



Bo 



B 4 - 



8Mv 
M 4 - s- 



-Be 



(18) 
(19) 

(20) 
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QRCS 



-Mt(I - I 1K )Im(B 1 + B 3 



AM 2 vI 1P 



Im(Si 
2Mt 



M 4 



B 3 ) 
lmB 6 



ImF, 



QRCS _ 



-M 2 t(I - I 1K )Tm(B 3 
4M 3 i'I 1 plm(B2 - B 4 ) 
2M 2 t(4:M 2 -i) 



B 4 



M 4 



-ImSfi 



(21) 



(22) 



We use next the forward kinematics and obtain for B n 
in the QRCS approximation: 



B QRCS 



M v Ff + tF 2 

{Mt{I - Ii K )Tm [Bl + B 3 + u(B 2 + B 4 )] 
+ 4M 5 W 1P Im [Bl - B 3 + i/(B a - B 4 )}} 
~t F x 



rF 2 



M 

+ 4M 2 vI 1P 4v 2 lm(A 3 



A 6 )} 



(23) 



where the invariant amplitudes of Lvov [Tsj are used in- 
stead of combinations of Bi appearing in the result. We 
make use of the well known physics content of the ampli- 
tudes A\ and A 3 + A 6 entering the final result of Ea. (|2*3^l . 
In terms of the hclicity amplitudes of real Compton scat- 
tering defined as T x > \> N ,\^\ N , 

M: 



W^ cs , one has 



Ax 



■A* = 



S ~ M2 T 1 i xl (24) 
MQVA/ 4 - su 2 ' 2 



2 x /Q 2 (s~M 2 ) 
1 



\v\jM 4 - su 



[Txl ji + T 1 i x : 

V 1 2' L 2 1 2' L/ . 



(25) 



It can be seen that the combination A 3 + Aq involves 
only helicity amplitudes without helicity flip, while Ax 
does flip photon helicity. Basing on the analysis of 
RCS |18| . we know furthermore that the real part of 
Re(^3 + A 6 ) is related to the forward nucleon polar- 
izability (a + f3), while KcAx is related to the back- 
ward polarizability [a — (3). It is interesting to observe 
that the backward RCS physics enters the expressions 
for forward ep-scattering amplitudes. This result is in 
fact model independent, since the contribution from the 
QRCS peak is always present in the full result. The role 
of the QRCS approximation used here is in neglection 
of the other contributions which might obscure this im- 
portant observation. Moreover, we see in Eq.JSBJ) that 
the forward combination of the RCS amplitudes enters 
B n multyplied by the integral Ixp which is peaked in the 
forward direction, as well, and the backward combina- 
tions with the backward peaked integrals Iq — Ixk- This 



result is in disagreement with the result of Ref. |14| . 
where the ep-scattering amplitude in the forward direc- 
tion is parametrized through the total photoabsorption 
cross section and the backward integral Iq. This discrep- 
ancy will be discussed in more detail in Section IVII 

In the following, we use the total photoabsorption cross 
section as the phenomenological input to Im(^43 + Aq) in 
the forward regime 18] 



lm(A 3 + A e )(s, t = 0) = ( - 



2v 



•so 



cip(0)-l 



(.26) 



with the total photoabsorption cross section a l °^ W 0.1 
mbarn. Furthermore, we adopted the standard Regge be- 
haviour in order to assure a correct energy dependence 
with the pomeron intercept ap(0) = 1.08 and the param- 
eter Sq — 1 GeV. For t dependence, we use an exponential 
suppression factor exp(Bt/2) with B — 8 GeV~ 2 , which 
provides a good description of Compton cross section for 
-t < 0.8 GeV 2 H3- 

The amplitude Ax in Regge regime is known to have 
the quantum numbers of a scalar isoscalar exchange in 
the i-channel which was successfully described by er- 
meson [Tsj or equivalently, by two pion exchange in the 
i-channel [l9|. Since the effective c-meson exchange 
does not result in a non-zero imaginary part in the s- 
channel, one should use the two pion exchange mech- 
anism accompanied by multiparticle intermediate state 
in the s-channel. In this work we estimate Ax through 
7T-ZV and pN contributions within the "minimal" Regge 
model for n and p photoproduction psl ] where reggeized 
description of high energy pion production is obtained by 
adding the t-channel meson exchange amplitude and (in 
the case of jp reaction) s-channel Born amplitude which 
is necessary to ensure gauge invariance. The reggeization 
procedure naturally leads to the substitution of each t- 
channel Feynman propagator by its Regge counterpart, 
V*(a*(t)), with 



i 

t—m'-i 



n = 



1 



sin 7ra,r (i) r(l + a w (t)) 



(27) 



with a v (t) — a4(i — m 2 ) and a'„ = 0.7 GeV -2 . Gauge 
invariance requires the s-channel piece to be reggeized in 
the same way, i.e., to be multiplied by (t—m 2 )T'^(a. R (t)). 
Here we list the results of the calculation of Ax'. 
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M(s - M 2 ) 



(28) 



m 2 p C 2 



M 2 



2M(s-M 2 ) {s-M 2 
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where CV 



C 



2y/2M 



2 v / 2Me 



u with f 2 NN /4ir = 0.08, and 
^ with fp^ = 0.103 0. Where 
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possible, and Q 2 were neglected as compared to s, M, 
m p in order to simplify the final expression. The scalar 
integrals Ao, Bq, Co for both irN and pN contributions 
are given in the Appendix. In Fig. [21 we show the ratio 
of the reggeized version of the integrals (A, B,C) P and 
the analytic results. The amplitude A± is defined as the 
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B n in the Regge regime 
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FIG. 4: Results for B n as function of Q 2 for different values 
of the LAB energy of the electron: 6 GeV (solid line), 9 GeV 
(dashed line), 12 GeV (dotted line) and 45 GeV (dashed- 
dotted line). 



FIG. 3: Ratio of the integrals with Regge propagators in the 
t-channel to the analytic results with pion Feynman propa- 
gators is shown for Aq (upper panel), Bq (upper panel), Cq 
(upper panel) as function of cm. scattering angle for beam 
energy E iab = 5.77 GeV. 

combination ^ [Bi+B 3 +u(B 2 +B4)], and the singularity 

j in the definition only cancels if taking the gauge invari- 
ant combination as described in (2^. The other feature 
of the result of Eas. H28l29|) is that both contributions are 

suppressed by factors jzfjps an d s _^ 2 , respectivelyly. In 
the case of the pion, it is interesting to observe this fact 
in view of somewhat surprising success of the effective 
field description of the SAMPLE data point on B n with- 
out pion contribution. This might give a hint that the 
pion contribution to B n at low energies is suppressed by 
the pion mass, if calculated to the same order. 



V. RESULTS 

We now present the results for B n using the model de- 
scribed in the previous section. In Fig. ^ beam normal 
spin asymmetry for the proton target is shown as func- 
tion of Q 2 for four different values of lab electron beam 
energy. In Fig. ^ the sum of forward A 3 + A 6 ~ a to t 
and backward A\ ~ 7T7r-exchange is shown. However, 
the impact of the latter is negligibly small in the kine- 
matics shown and is far below 1% for all energies consid- 
ered. The reason for that are the both suppression factors 
(s_M^) an d ( s ?M-i) together with the Regge suppression 
due to the pion trajectory, as compared to the Pomeron. 



So, the log 2 m enhancement of B n is absolutely irrele- 
vant in this extreme forward kinematics. Recalling that 
this log 2 m term factorizes the backward Compton am- 
plitude contribution, this is by no means a surprise, since 
it is unphysical for bacward mechanisms to dominate for- 
ward observables. This result suggests in turn that the 
QRCS approximation should not be expected to work in 
the forward kinematics. In fact, this effect was observed 
on the example of irN intermediate state contribution in 
[Icj . However, if going to backward angles, this contri- 
bution increases as ~ \t\ ln 2 (|i|/m 2 ) and, in general, the 
approximation is justified. In the forward kinematics, 
the energy dependence in Ea. (|23[) originates thus mostly 
from the term ln-^=, so that the asymmetry increases 

logarithmically with the energy. This is an interesting 
observation and is on contrary with the dominant en- 
ergy dependence of B n at low energies where it decreases 
with energy as m/E. It furthermore gives good outlook 
for measuring B n through high energy region at forward 
angles. Still, a full calculation (beyond the QRCS ap- 
proximation) is needed to confirm this behaviour. Com- 
paring to the calculation of |14l |. the quite different way 
the real Compton scattering amplitudes enter the final 
result for B n leads to inversed energy dependence and 
consequently, the largest disagreement between the two 
calculations amounts in an order of magnitude for the 
JLab energy, lowest among those considered. This dis- 
crepancy will be discussed in details in the next section. 
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VI. DISCUSSION 

As it has been noticed before, the presented calcula- 
tion is in disagreement with the result of Ref. (TiJ ■ Their 
approach consisted in taking only the part of doubly vir- 
tual Compton tensor which survives in the exact forward 
limit of real Compton scattering (see Eq. (14) of [T^l 



{-(PK) 2 g^- 
+ {PK) [P»<fi 



(<?i? 2 )W 



(29) 



with A the forward amplitude of doubly virtual Comp- 
ton scattering which is then related to the proton struc- 
ture function F\ and the total photoabsorption cross sec- 
tion <t*°*. The forward spin independent doubly virtual 
Compton scattering tensor is usually written in terms of 
the DIS structure functions W\, W 2 , 



or 



1 

~Mv 



(?1 • 92) 

Pen 



W 2 {v,i 
Wi 



2 ) 



(91 ■ 12 



Lt 



pv 



Pq 2 
(91 • Q2)' 



(30) 



Qi W 2 , 



where the tensors should be taken in the limit qi — 
q 2 . In the scaling limit, one has W\(y,Q 2 ) — > F\{x), 
W 2 (u,Q 2 ) — > F 2 (x), and the Callan-Gross relation im- 
plies F 2 = 2xF\ . We rewrite now the tensor of Eq. (f 3()|) 
in terms of the DIS tensors Oy V , O^- 



= [PKf {or - A 



(31) 



We proceed by projecting this tensor onto the basis of 

Prange, make use of ^^f 2 - = 2 + ^g? 2 , and obtain for 

the amplitudes' combination which multiplies the In 2 m 
term: 



B x + B 2 + v(B 2 + B 4 



{PK) 2 Q 2 + Q 



2M 2K 2 



A, (32) 



The kinematical limit of Ref. [l4j corresponds to ne- 
glecting t comparing to the photons' virtualities. Then, 
qi = q 2 = K, and the right hand side is non-zero. In- 
stead, throughout this work, the non-forward RCS kine- 
matics with the on-shell photons is used, and the contri- 
bution of the cross section (i.e., helicity conserving helic- 
ity amplitudes) to Eq. l|32|l vanishes. It is an interesting 
phenomenon, that the result of taking the limit Q 2 2 — > 
and i — s- is order depending. 

The dependence on the order of taking the limit Q 2 — > 
and low photon energy limit has been observed before 
in the context of polarizabilities of RCS and generalized 
polarizabilities of VCS [24|, and even for the leading term 
in photon energy of the forward RCS vs. forward doubly 
VCS amplitude F(v,Q 2 ) 



lim lim F(v,Q 2 ) = 
lim MmFlv.Q 2 ) = 

Q3_ > 0f->0 



M 



■e N (e •£), 



(33) 



M 



n N (2e N + kn)(s* ■ e), 



with ejv and kn the charge and the a.m.m. of the nucleon 
(1 and 1.79 for the proton and and —1.91 for the neu- 
tron, respectively). So, it is not completely unexpected 
that the two different limits do not match. In Fig. [31 
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FIG. 5: Mandelstam plot for doubly virtual Compton scat- 
tering with Q\ = 0.3 GeV 2 and Q\ = 0.01 GeV 2 . For further 
details, see text. 

the Mandelstam plot for the general case of doubly vir- 
tual Compton scattering is shown for Q\ = 0.3 GeV 2 
and Q 2 , = 0.01 GeV 2 . The plot shows different kinemat- 
ical regions on the plane v = 52± and t. The energy 
threshold for the reaction in the s-channel, s = M 2 , is 

represented by the straight line v = — ^ M 2 , so that the 
s-channel reaction (direct Compton scattering) is possi- 
ble to the right of this line. The requirement that the 
scattering angle takes physical values, restricts allowed 
values for the variables v and t to the area between the 
two curves depicted in the right half-plane. The upper 
line corresponds to cos = 1 (forward scattering) , while 
the lower one ot cosO = — 1 (backward scattering). The 
physical scattering can only occur inside this region. The 
forward line crosses the t axis at t = —{\J~Q\ — \/Qf) 2 , 
while the backward line at t = — (-\/Qi + \/ r Q 2 ) 2 - The u- 
channel (crossed Compton reaction) kinematical regions 
are obtaine by v — > — v. The variable v should not be 
confused with the energy of the virtual photon, which is 
sometimes denoted as v, as well. In the lab frame, the 

initial photon energy is given by lo\ = — — 2 I M +< ^ 1 1 while 
" — — — — — , with w = (v + oi ) the m- 



4M 



variant mass squared of the intermediate hadronic state. 
For each value of the external variables s and t which 
one measures in a ep-scattering experiment, the inte- 
gral over the intermediate electron phase space involves 

$M2 dw 2 and J dQ\dQ\. The latter integral is per- 
formed over the ellyptic areas shown in Fig. [21 unambigu- 
ously defined by s, t and w 2 . In terms of the Mandelstam 
plot in Fig. [SJ this corresponds to integrating over a se- 
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ries of such Mandelstam plots for each pair of Q\ , Q\ at 
the value of t fixed by the experimental kinematics. In 
the special case of Q\ = Q 2 which is always contained 
within the integral, as can be seen in Fig. |2 the Mandel- 
stam plot is deformed, so that the forward line becomes 
t = 0. If Q\ 7^ Q 2 , however, t — is only approached 
asymptotically for v — ► 00. This is the reason why the 
limit t — > is not well defined for the general case of 
the doubly virtual Compton scattering: this limit takes 
one outside of the physical region of the reaction. If one 
preforms the limit Q\ 2 — > first, keeping t 7^ 0, one en- 
ssures that the line t = is contained within the physical 
region, and the limit can be performed. 

Technically, the difference between the two limits orig- 
inates from treating as small quantity either the momen- 
tum transfer (forward doubly VCS) or the virtualities of 
the photons (RCS). We note that the squared logarithm 
enhancement occurs if the hadronic amplitudes is non- 
zero in the exact limit of hard collinear photons. There- 
fore, the leading (in log 2 m) contribution should corre- 
spond to Compton scattering with the photon virtualities 
of the order Q\ = — (k — ki) 2 w 2mE and the momen- 
tum transfer of the ep reaction which one measures on the 
experiment. It means that in order to justify the limit 
of one nas to go to external momentum transfers 
smaller than 2mE which is of the order of 10 -3 GeV 2 in 
the kinematics considered here. Apart from difficulty in 
achieving this kinematics experimentally, it is suppressed 
by the kinematical factors in the general expression of 
Eq.JHJfor B n , which ensures this observable to vanish in 
the exact forward direction. Basing on these considera- 
tions, we see that in order to obtain adequate estimate 
of the double log effect for B n , one has to split the full 
doubly virtual Compton amplitude into the non-forward 
real Compton amplitude which gives the desirable effect, 
plus the rest which is regular in the QRCS point and can 
at most give a In to effect. 



VII. SUMMARY 

In summary, we presented a calculation of the beam 
normal spin asymmetry in the kinematical regime of high 
energies (i?; a fc = 6 — 45 GeV) and low momentum transfer 
to the target, Q 2 < 0.5 (GeV/c) 2 . This observable ob- 
tains its leading contribution from the imaginary part of 
the two photon exchange graph times the Born amplitude 
and is directly related to the imaginary part of doubly 
virtual Compton scattering The resulting loop integral 
obtains a large contribution from the kinematics when 
both exchanged photons are nearly real and collinear to 
the external electrons. We adopt the QRCS or equivalent 
photons approximation which allows to take the hadronic 
part in the external kinematics out of the integral and 
to perform the integration over the electron phase space 
analytically. For the hadronic part, we use the full real 
Compton scattering amplitude and show, that both for- 
ward (no proton and photon helicity flip) and backward 



(photon helicity flipping) amplitudes do contribute in this 
observable in the forward kinematics. For the forward 
Compton amplitude, we use optical theorem as the in- 
put. For backward Compton amplitude, we provide an 
cxplicite calculation which is due to two pion exchange 
in the t-channel. The resulting values of the asymmetry 
are in the range 4 — 8 ppm for the energies in the range 
6 — 45 GeV. We find furthermore that the double loga- 
rithmic enhancement does not dominate B n in forward 
regime since it comes with helicity-flip Compton ampli- 
tude which highly suppresses this behaviour. Finally, we 
make an interesting observation that the non-analyticity 
of the doubly virtual Compton scattering leads to the 
order dependence of taking limits t — > and Q\ 2 ~ * 0> 
which is a generalization of the known problem of match- 
ing the low energy limits of real and doubly virtual Comp- 
ton scattering. 
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VIII. APPENDIX A: INTEGRALS OVER 
ELECTRON PHASE SPACE 

In this section we present calculation of the integrals 
over electron phase space appearing in the QRCS approx- 
imation. First we calculate the scalar integral Iq, 



h = 



k\dk\ 



o 2£i(2tt 3 ) J (fc-fci) 2 (fc'-/ci) 2 ' 



(34) 



where the upper integration limit kthr corresponds to 
the inelastic threshold (i.e. pion production), k t hr = 
^J (s-(M+m p i) ) _ ~ 2 m nex t introduce integration 

over the Feynman parameter = r^q^r^p ■ We 
chose the polar axis such as k\-(k—xq) = k\\k—xq\ cos 0i 
with \k — xq\ 2 = k 2 + x(l — x)t and perform angular inte- 
gration. We furthermore change integration over dk\ to 
integration over dimensionsless z = E\/E 



-8ttH 



dz 



e 2 t y 1 Z > 



(35) 



•In 



m 2 4m 2 

E 2 t 



(l-z) 2 + Jz 



E 2 



m 2 4m- 

E 2 t 



(l-z) 2 -Jz 2 - 



E 2 



To perform the integration over the electron energy, 
we follow here the main details of the calculation in Ap- 
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pendix A of Ref.Q. The result reads 
1 



h = 



-32tt 2 



' i 2 E 2 



where Sp(x) is the Spence or dilog function, Sp(x) — 
- § ln(l - xt) with Sp(l) = tt 2 /6. In the high energy 



limit. 



1, we recover the result of Ref.[lJ|. 



We next turn to the vector integral 

^ kldkx f k>i 



Jf = 



o 2^(2^) 
hpP* 1 + I 1K K» 



1 (fc-fci) 2 (fc'-/ci) 2 



(37) 



It cannot depend on due to the symmetry of I\ 
under interchanging k and kl . To determine these two 
coefficients we have a system of equations, 



k°, 1k + i-a, - if - JL, / £ 



(38) 



(27r3)S 1 Q 2 



Using the same approach as for Iq, we obtain for 7° 

^ - E thr 

E f E zdz 



II = 



-8tt 2 £ 



,2 m 2 4m^/i _ \ 

£ 2 t I.J- z; 



(39) 



■In 



2 m 2 4m 



f(i-,) 2 + j- 



,2 _ m^ 
"E 2 " 



2 m 2 4m 



-47T 2 t £ 



In 



+ 1 



V—t E t hr 

m E 
Eth 



E 



In 1 



Ethr 

E 



E 2 

(40) 
(41) 



Finally, we consider the integral 1% where we perform 
angular integration, 



h = 



1 



d 3 k x 

(2^)y -E~Ql 



(42) 



dz In 



m_ , ;i m 2 , /„2 ™1 

2 B 2 V E 2 \l ~ E 2 



m- /-1 m 2 / 2 _ m. 



After integrating by parts and changing variables con- 
sequently z = ^7 cosh y and j/ = hit, we arrive to 



-Ethr , 2E t hr 



U = 27r^ln 



2tt 1 



Ethr 



(43) 



In 1 - 



Ethr 

E 



In these last two integrals, it is important to keep -^g^ 
unequal to 1 till the end to ensure the convergence of the 



integral. Solving the system of linear equations for the 
coefficients, we obtain: 

8 - M 2 E t hr , 2E 
hp — TT-r-n ln- 



M 4 - su E 

1 APK 
Ixk - Q>h -QThp. 



(44) 



IX. APPENDIX B: SCALAR INTEGRALS FOR 
HELICITY FLIP AMPLITUDE 

The vector and tensor integrals can be reduced to the 
scalar ones by means of standard methods [25J. The re- 
maining integrals to be calculated are the two, three, and 
four-point scalar integrals. Here we are only interested 
in the imaginary part of these, therefore there are only 
three integrals with non-zero imaginary part: the two- 
point integral 



OS 



Im 



1 



1 



(2tt) 4 pi - ml {P + K - p^ 2 - M 2 



8tt Js 



(45) 



with Ip^I 



_ /(s-M 2 +m 2 ) 2 



m 2 ; the three-point one 



d^Pir 1 



1J " : 1 (2n) 4 pi - ml (k - p^) 2 - ml 



{P + K -Pn) 2 -M 2 



•In- 



2E„ 



8tt(s - M 2 ) 
and finally, the four-point integral 
d A Pn 1 



(46) 



A n = Tm 

1 (27r) 4 (fc-p.) 2 -m 2 p 2 -ml 



1 



1 



(fc' - pn) 2 -ml(P + K- p n ) 2 - M 2 
1 



8vrQ 2 (s - M 2 +m 2 ) 
1 



-. 4m 2 E 2 
1 + Q 2 P 2 



hi 



1 + ^ + 1 



-. 4m' 2 E 2 -. 

t H 1 



(47) 



These integrals should however be reggeised as de- 
scribed in Section by substituting the Regge propaga- 
tor instead of the Feynman one. Denoting t\ — (k — p n ) 2 
and ^2 = {k' — Prr) 2 the momentum transferred by the 
pions in the i-channel, we have for the reggeized version 
of scalar integrals: 



(CS) 



1 Pw 

32tt 2 



dQAh-mi^ia^h)) 
■{t 2 -ml)V*{a„{t 2 )) (48) 
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= -J—^ / dn^h-ml) (49) M = Im f - 



p£K(ti))p*K(t2)), i P i 



(fc' - p p ) 2 - m 2 (P + K - Vp ? ~ M 2 
1 1 



{A ° )R = 3^^/^^K(*0)^K(*2)) (50) = g^Q^M^g) ^1 + ^ 

Similarly, in the case of the p-exchange in the s- ( f\~> ^ 

grals with non-zero imagii 

d^Pp 1 1 



i Q 2 p 2 

channel, the integrals with non-zero imaginary part are: . v p j (51) 



C£ = Im 



1+7^,-1 



(2tt) 4 p 2 p - m 2 p (P + K - p p ) 2 - M 2 

UpA 

8tt Js' 



Im 



Pp 



(27r) 4 JJ 2 _ m 2 (fc _ p p )2 _ m 2 



(P + K - p p ) 2 - M 2 



8tt(s - M 2 



In 



2£ p (s - M 2 
Mm 2 



with |pp 

££ p (m 2 - mi 
Mm 

mass, a 



(s-M 2 +rn 2 p ) 2 
Is 



i 2 and a 2 



E 2 m 2 p 



2s 



"+" i( m p — m V) 2 ■ If neglecting the pion 
The reggeization procedure is the same 



as for irN- intermediate state. 
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